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Abstract 

We show how generalised unitarity can be used to determine the one-loop dilatation oper¬ 
ator in J\f = 4 super Yang-Mills. Our analysis focuses on two sectors, namely the bosonic 
SO(fo) sector and the SU( 213) sector. The calculation is performed on shell, with no off- 
shell information introduced at any stage. In this way, we establish a direct connection 
between scattering amplitudes and the dilatation operator of the A^ = 4 theory. 
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1 Introduction 


In this paper we continue the study initiated in [1], whose ultimate goal is to relate scat¬ 
tering amplitudes in Af = 4 super Yang-Mills (SYM) theory, and their properties, to the 
dilatation operator in the same theory. In a sense, we retrace the history of the develop¬ 
ments in the calculation of loop amplitudes triggered by Witten’s twistor string theory [2] : 
In [1] we employed MHV diagrams [3] at loop level [4] in order to compute the dilatation 
operator at one loop; here, we proceed to apply a powerful variant of unitarity [5,6] known 
as generalised unitarity [7,8] which, as we shall see, allows for an even more efficient cal¬ 
culation of the dilatation operator. As generalised unitarity turned out be more practical 
than loop MHV diagrams, we will see how our use of generalised unitarity will further 
simplify the already remarkably simple calculation of the dilatation operator performed 
with MHV rules. 


The use of unitarity in deriving the dilatation operator is welcome also from a conceptual 
point of view, since the only ingredients of the calculation are on-shell amplitudes - with no 
off-shell information being introduced. This supports the hope that using this approach one 
may be able to connect directly the amplitudes and their hidden structures and symmetries 
to the integrability of the dilatation operator in Af — 4 SYM. 

At one loop, the no-triangle property [5] of the one-loop S-matrix of Af = 4 SYM implies 
that maximal cuts employed in [8] are enough to completely determine all amplitudes of 
the theory. Similarly, we identify certain quadruple cuts which are sufficient to determine 
the dilatation operator at one loop. The reason why this is correct lies in the simplicity 
of the object under consideration, namely a two-point function, or a single-scale object 
in momentum space. At one loop there are precisely four fields to be connected, which 
explains why quadruple cuts are enough. 

In more detail, we will focus on the dilatation operator in the 50(6) and 577(213) sectors 
of Af = 4 SYM, which we will derive by computing the two-point functions (0(xi)0(x 2 )) 
of the appropriate composite operators. The 50(6) sector was studied in [9], where the 
connection to integrable spin chains was first made, and is closed (only) at one loop. The 
50(213) sector is closed to all loops and was considered first in [10]. At one loop and 
in the planar limit, only contractions of pairs of adjacent fields survive. In all cases we 
consider, the result of our calculation turns out to be proportional to the quadruple cut 
of the same integral appearing in the original one-loop calculation performed by Minahan 
and Zarembo. This integral is given by 


I(x 12 ) — Jd D z A 2 (aq - z) A 2 (x- 2 - z ), 

where x i2 := X\ — x 2 and 

A(x) := 
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is the scalar propagator in D dimensions. In momentum space, it appears as the Fourier 
transform of the simplest single-scale integral, namely a double bubble, 
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Figure 1 : The double-bubble integral relevant for the computation of I(x 12). 

where L := L± + 1 / 2 - By using the quadruple cut where the propagators with momenta 
L 1 , L2 '■= L — L\, L% and L 4 := L + L 3 are put on shell, we will be able to identify the 
coefficient of this double bubble in all relevant cases - without ever performing an integral. 
The cut double bubble can then be lifted to a full integral, and by picking its ultraviolet 
(UV) divergence I(x 12)| uv , 

/(a; i 2 )l UV = 7' 8 ^ ' (4n 2 x 2 12 ) 2 ’ 

we can immediately write down the dilatation operator. Let us also mention that other 
applications of unitarity to the calculation of n-point correlators and correlation functions 
of Wilson lines have appeared in [11-13]. It is worth stressing the two key reasons why 
generalised unitarity is particularly powerful for the case considered here, namely that of the 
two-point function. First, as we have already mentioned, quadruple cuts are precisely the 
right set of cuts to identify the relevant loop integral; furthermore, the relevant integrals are 
guaranteed to have a single scale. The cut integral can then be lifted to a full loop integral 
without introducing spurious discontinuities, in complete analogy to the case of splitting 
amplitudes studied in [14], Finally, we also note that our approach to the computation of 
the dilatation operator differs from that of [15,16] in that no infrared divergences appear 
at any stage in our calculation. 

The rest of the paper is organised as follows. In the next section we use generalised 
unitarity to obtain the dilatation operator in the SO ( 6 ) sector. In Section 3 we move on 
to the SU(2\3) sector. This case is particularly interesting as it involves fermions as well 
as scalars. There are several contributions to consider and the structure of the dilatation 
operator is more elaborate than in the pure scalar sector, hence our tests are more stringent. 
Finally in Section 4 we compare recent on-shell and twistorial approaches for the calculation 
of the dilatation operator, and also make a few suggestions for future work. 


2 The dilatation operator in the 50(6) sector 

In this section we will compute the dilatation operator of J\f = 4 SYM in the 50(6) sector 
using generalised unitarity. This calculation was recently performed in [17] and [1] using 
MHV diagrams in twistor space and momentum space, respectively. Here we depart from 
these off-shell approaches in favour of a fully on-shell calculation. 

Operators in the 50(6) sector have the form 

O AiBi,A2B2,...,AlBl (0 ■ r Tf’(0A 1 S 1 Od ' ' ' ^AlBl (^)) ■ (^-l) 
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At one loop and in the planar limit, it is sufficient to consider contractions of pairs of 
adjacent fields (in colour space). The relevant part of each operator is then 


0 (x i) = • • • ct> a A B (xi)<t> b c D (x r) • • • {T a T b r 3 , 

0(x 2 ) = ■■■ d> C A'B'(x2)ct> d c, D '(x 2 ) • • ■ ( T c T d y m . 

The calculation is then effectively equivalent to that of the following two-point function 

((0AB0CT))( x l)(0A , S'0C'D')( a; 2)) , (2.3) 


whose expected structure is 

(( 4>AB4>CD)(xi)((f>A'B , (f>C l D l )(x2 )) = ■A^ABCD^A’B’CD’ + BeABA’B'tCDC’D' +CtABC'D'^A’B’CD ■ 

(2.4) 

These three terms are usually referred to as trace, permutation and identity. The dilatation 
operator can then be read off from the UV divergences of (2.4), hence we only need to 
compute the UV-divergent parts *4.uv; £>uv ; Cuv of the coefficients A , B and C. These are 
expected to be equal to [9] 

Aw = -, Buy = —1, Cuv = 1 ■ (2.5) 

As in [1], we choose the following SU(A) assignments in (2.2) as representatives of these 
three flavour structures: 



ABCD 

A’B’CD’ 

Tr 

1234 

2413 

P 

1213 

3424 

1 

1213 

2434 


For each case there is a single cut diagram to consider. The integrand is constructed with 
four cut scalar propagators with momenta L i: % — 1,... ,4, and one on-shell amplitude, as 
shown in Figure 2. The operators are connected to the amplitude via appropriate form 
factors, which in the scalar case are simply 


:= fd l xe iL * 

= (2tt) 4 A C4) (L-h- i 2 ) 5 aa '6 w , 

where we have used 0 and <fi to denote two scalar fields having distinct A-symmetry indices 
as is sufficient for our purposes, see (2.6). Note that the ii represent the on-shell (cut) 
versions of the loop momenta Lj. 

The relevant amplitudes for the three flavour assignments considered in (2.6) are: 1 


Tr : 

A(l^ 12 4 9il3 3^ 24 2^ 34 ) = 

(13}<24> 

(2.8) 

(12)(34) ’ 

P : 

A^ 42 ,4^4, S^ 34 , 2 ^ 3 ) = 

-1, 

(2.9) 


1 In the following expressions we omit a factor of 3 ym> which 
calculation. 


will be reintroduced at the end of the 
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Figure 2: The single cut diagram contributing to the dilatation operator at one loop. 


1: A( i^w^ 34 ^ 4 ^ 13 ) 


(13)(24) 
(23)(14) ' 


( 2 . 10 ) 


Three observations are in order here. First, we note that the same integrands as in the 
approach of [1] have appeared, with the important difference that, in that paper, the spinors 
associated with the on-shell momenta are given by the appropriate off-shell continuation 
for MHV diagrams. Here the spinors for the cut loop momenta do not need any off-shell 
continuation. Furthermore, for the case of the P integrand there is obviously no difference 
between the two approaches, and the resulting integral is given by a double bubble where 
all the four propagators are cut. In the other two cases, this integral is dressed by the 
appropriate amplitude. Finally, we note that the colour factor associated with all diagrams 
is obtained from the contraction 


... (T 6 T“)h • • ■ Tr (T a T h T c T d ) ■ ■ • ( T d T c ) l m ■■■ = ■■■ N 2 8^5) • • • , (2.11) 


where the trace arises from the amplitude and the factors • • • (T b T a ) l j ■ ■ ■ and • • • ( T d T c ) l m ■ ■ ■ 
from the operators (and we indicate only generators corresponding to the fields being con¬ 
tracted). We now proceed to construct the relevant integrands. 


The trace integrand 

In this case the relevant amplitude (which multiplies four cut propagators) can be rewritten 


/1Q\ / r )A\ TV (0.0-0.0.\ 0(0. 



( 2 . 12 ) 


where we have used l\ + £2 = — (^3 + ^ 4 ) : = L. Having rewritten the amplitude in terms 
of products of momenta, we lift the four cut momenta off shell. The resulting integral has 
the structure of a product of two linear bubbles, 


2 r d D u L 1 ; rd D L 3 l 3)1 

L 2 ./ (2n) D L\ (L - L x ) 2 ./ (2vr)^ L\ (L + L 3 ) 2 ' 


/ 


/ 


(2.13) 


Using the fact that 



(2.14) 


2 We define Tr + (abed) := ( ab) [be] (cd) [da]. 
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where (in Euclidean signature) 


Bub(L 2 ) 


f d D K i i r(2~ f)r 2 (| -1 ) 2 g _ 2 

J (2n)° K 2 (K + L) 2 ( 4 vr)f T(£> - 2 ) 1 


(2.15) 


we find that (2.13) is equal to 1/2 times a double bubble. Using (1.4) we finally get 
*4.uv = 1/2. Note that in arriving at this result we have performed a Fourier transform to 
position space using 


D ip-x 


d D p 
{2n) D (p 2 ) s 


r(f-s) i 

4 s 7r(s) (x 2 )^~ s 


(2.16) 


Furthermore, in the definitions of Aluv, ^uv, and Cuv (and in general of any other UV- 
divergent coefficient in the rest of the paper), a factor of A/( 87 r 2 ) x (l/( 47 r 2 a; 2 2 )) x (1/e) 
will always be understood, with A := gy M N. 


The P integrand 

No calculation is needed in this case, and the result is simply given by minus a cut double¬ 
bubble integral. Lifting the cut integral to a full loop integral we get B\jy — — 1. 


The 1 integrand 


The relevant amplitude in this case is 

(13)(24) 
(23)(14) 


(12)(34) 
(23)(14) ' 


(2.17) 


Thus the first term in (2.17) gives the cut double-bubble integral, whereas we can use 
on-shell identities to rewrite the second term as 

(12)(34) (12)(34)[34] L 2 

(23) (14) (23) (14) [34] 2(«, • 64 ) ' y ' ' 

Lifting the cut propagators of the second integral to full propagators, it is immediate to 
see that this term produces the integral represented in Figure 3. This integral is finite in 
four dimensions and thus does not contribute to Cuv- We then conclude that Cuv = 1. 



Figure 3: The finite integral corresponding to the term in (2.18). This integral is irrelevant for 
the calculation of the dilatation operator. 


A comment is in order here. In principle an ambiguity is still present corresponding to an 
integral such as that of Figure 3 but with one of the four propagators L 1 ,..., L 4 collapsed 
(say Lfi), which is UV divergent. This integral can be excluded by looking at a triple 
cut corresponding to cutting the propagators Lj, L 3 as well as the middle propagator in 
Figure 3. 
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For later convenience, we explicitly write down the form of the UV-divergent part of 
the correlator (2.4), 


(( 4>AB(pCD ) (^1 ){<j>A'B' <t>C'D') l m ( x 2)) 

J uv 

= ~ ' ^2 (^ 2 ( 2 ' 1 2)^m^) (2 e ABCD^A'B l C l D' ~ C ABA'B'^CDC'D' + ^ABC'D'^A'B'CD^j 

In terms of a spin-chain Hamiltonian, this can be represented as [9] 

A /I 

H 2 = Tr + 1 — P 

8tt 2 V2 


(2.19) 


( 2 . 20 ) 


3 The dilatation operator in the SU( 2 | 3 ) sector 

In this section we consider the closed SU( 2|3) sector. This is particularly interesting, as it 
involves also fermions. Indeed, operators in this sector are formed with letters taken from 
the set {^i Q , 0m}, with a = 1,2 and A = 2,3,4. We thus have one fermion and three 
scalar fields. The dilatation operator in this sector was derived in [10]. Its expression is 
given by 



(3.1) 


where in this notation the three scalar fields are labelled by A, B = 2,3,4. In the following 
we are going to rederive (3.1) using an application of generalised unitarity. 

As for the SO (6) case, in the planar limit only contractions between nearest-neighbour 
fields in 0(x 1 ) and 0(x 2 ) have to be considered. The first two terms on the right-hand side 
of (3.1) denote the scalar identity 1 and permutation P structures already familiar from 
the SO( 6) case (the trace structure is absent given the restricted choice of scalar letters). 
The novelty is that now we have to consider two additional types of contractions: scalar- 
fermion —>■ scalar-fermion, and two-fermion —>■ two-fermion, as indicated in the remaining 
terms in (3.1). 


Scalar-fermion —>■ scalar-fermion 


In this case we are interested in a fermion field and one of the scalars <fi 12 , </>i 3 , or 0i4. 
Without loss of generality we will consider 0 12 . There are two cases to consider, 


U : ((0i2'0io)( a; i)('^234 0 :^ 34 )(^ 2 )) j (3-2) 


and 

S : (034^2346)0*2)) , (3-3) 

where the letters U and S indicate whether the contractions between the two fields are 
unswapped or swapped. The relevant form factor is 




MV 


Ac 

l 2 


;£) := 


e iL x <0|(«‘„)(i)|^(< 1 ),<(< 2 )> 

(27r) 4 ^ (L-O- i 2 ) Xl 8 aa '5 bb ', 


(3.4) 
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and similarly for 0 (x 2 ). 

We begin by considering the U case. By contracting the two form factors with the four 
planar permutations of the full amplitude, we obtain'* 


x 


A 2 A 3 5 aa' 6 bb' 8 cc' 5 dd' 

2^ 1 4^ 34 ) Tr( T a 'T b ' r j lc 'T d ''\ _|_ 12 2 ,| b 4^34 3^234 ^ Xr(T a, T b, T d, T c 


_ S^ 234 4^ 34 2^’ 1 ) r Yx(T a 'T c 'T d 'T^' _ j4(1^ 12 4^ 34 3^ 234 2 1 A)Tr(T ,a, T d 'T c ‘T b ') 


(3.5) 


At large iV there is only one leading contribution, corresponding to the term with the 
amplitude 


A(l^ 12 4^ 34 3^234 2^!) 


(13) (34) 

(14) (23) ' 


(3.6) 


It is given by 


-^(l^ 12 ,4^ 34 , S^ 234 , 2^) \ 2 p\* 


(4 £\ h)pp 
2{£i • Q 


N, 


pp- 


(3.7) 


The cut integral to consider is thus 

hi,-fd l i 1 d a i 3 s <+ \il)s^\tl)sM([L-hf)s^((L + i 3 f) ■ n^, ( 3 . 8 ) 

where by Lorentz invariance must have the form 

Ipp = ALpp . (3.9) 

A simple PV reduction shows that the UV-divergent part of the coefficient A is equal to 3 4 
Auv = 1/2. 

For the S case, we get the single leading contribution to be 

-A(iA 2 ,4^234,3<>34 )2 yi )A 2 A 4 = = ^ (3. 10 ) 


The relevant integral is now 

4 s := fd%d%S^((i) 5 i+ \ll) 5 ^((L-l,) 2 ) <S<+> ((L + < 3 ) 2 ) • N 0i) 

= ALpp , (3-11) 

where a PV reduction shows that A = —1/2. Note that in arriving at this result we have 
discarded finite integrals, which do not contribute to the anomalous dimensions (more 
precisely, in all calculations the only other finite integral appearing is the kite, depicted in 
Figure 3). 

3 Two out of the six possible contractions obviously do not contribute at large N. 

4 We recall that we omit a factor of A/(87 t 2 ) x (1/ (47r 2 a:f 2 )) x (1/e) after Fourier transforming to position 
space. 
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Summarising, the scalar-fermion —y scalar-fermion case gives ±1/2 times a double¬ 
bubble integral, for the U/S case, respectively. This has to be compared to the tree-level 
expression 


f d D Li L\pfi 
J (2n) D Lf(L — Li) 2 


1 

2 


%Bub (L 2 ) . 


(3.12) 


Thus for the two-scalar two-fermion case we get: 


1:1, P : -1 , 

and the corresponding contribution to the spin-chain Hamiltonian is 5 



in agreement with the corresponding terms in (3.1). 


(3.13) 


(3.14) 


Two-fermion —> two-fermion 

In this case we consider the four-point correlator 

(Wa^l ■ 

The form factors of 0(x i) are given by 

J d ‘ xe ‘ Ll <olW±s)(^OI«'(fi).«f2)> 

= (2n)V a \L-t I -e 2 ) -i(Api- AjA 2 J ab ’6 ba '), 


(3.15) 


(3.16) 


and similarly for the form factor of 0(x 2 ). Note the factor of 1/2 appearing because of the 
presence of two identical particles in the state. Contracting the two form factors with the 
four planar permutations of the full amplitude, we get 


--(A* A l5 aa '5 w - AjA 2 J ab 'S ba ') (Xl\*5 cc '5 dd ' - X^X^8 cd ' 8 dc ') 
x A{Y ^ 2^ 1 3^ 234 4 ^’ 234 ) Tr(T a/ T b ' t c 'T d '\ _ AiV^ 1 2^ 1 4^’ 234 3 1 /’ 234 ) Tr^^T^T^T^) 
±r 4 (l^’ 1 3^ 234 4^ 234 2 ^ > 1 ) r £x(T a 'T c 'T d 'T b ') — A^Y ^ 1 4^234 2^ 4 )Tr(j ia/ 

(3.17) 

In the large -N limit we only need to keep the following terms out of those in (3.17): 


AL(r /,1 ,2 v ’ 1 ,3 ,/ ’ 234 , 4 te4 ) AJA^A* ± 2^, T^ 234 ,3 V ’ 234 ) AjA^A|A 


A(V h , 3 V ’ 234 ,4 fe 4, 2^ )X 1 a X 2 p X 3 $ Xi - H(l^ 1 ,4^ 234 , 3^ 234 , 2^) A* A^A^A 


(3.18) 


5 Here we also reinstate powers of gy M from the tree-level amplitudes, of N, arising from colour con¬ 
tractions, and a factor of l/(87r 2 ) arising from the UV singularity (1.4) of the double-bubble integral 
(1.3). 
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where the relevant four-fermion amplitudes are 


2^! 3h234 ^V’234^ _ 

A(V>h 2^ 4^234 3^234 ^ _ _ 

j 4 ( 1 ^’ 1 3^ 234 4^234 2hi j = 

yl(l ^’ 1 4 ^ 234 3^234 2 ^^ — 
Using (3.19), we can rewrite (3.18) as 


<34) : 


(23}(41} ’ 


(34)' 


(24)(31) ’ 


(34)' 


(13)(42)’ 
(34) 2 


1 

4 


(14)(32) ' 

+ (^ 1 ^ 2 )q/ 3 (^ 3 ^ 4 )q^ 


2(fe • 4 


+ (-i -H- 


(3.19) 


(3.20) 


The term with td £2 is simply a relabelling of the integration variables, and we conclude 
that the one-loop integrand is given by 


(hh)a/3 (44)q/ 3 + (44)a/3(44) a$ 


2(^2 • h) 

Thus we have to consider the cut-integral 


= N. 


a(3al3 


(3.21) 


L a/36i(3 


U%d%5^ + \e i )5^\^)5^{{L-^) 2 ) 5^{{L + ^f) ■ N af}&$ . (3.22) 


It depends on only one scale L , hence it has the form 

4/36/3 = ^ -^ 2e «/? e a/3 + B (LaaLpp + L^Lp^) 


(3.23) 


Contracting (3.22) and (3.23) with e Q/3 e°^ and ( L aa L^ + L^ a L al3 ) we can solve for the 
coefficients A and B. The result for the corresponding UV-divergent parts is 


Ajv — 0 , Buy — 1/6 . 


(3.24) 


At this point we lift the four cut propagators to full propagators, so that the cut double 
bubble becomes a full double-bubble integral. The conclusion is then that the UV-divergent 
part of the integral representing the two-fermion —y two-fermion process is a double bubble 
with coefficient 


g (LaaLpp + B a$Lpci) ■ 

This result has to be compared with the planar contractions at tree level, 

.tree = f * D I* L la$ (L - 

<*(>*& ' ] (2v t) d LKL-Li ) 2 ' 


(3.25) 


(3.26) 


After a similar PV reduction of the L\ integration in (3.26), we find that is given by 

a scalar (single) bubble with coefficient 


1 

6 


( - L 2 e aP e 


a(3 


+ L a$ L Po‘) 


(3.27) 


9 













which using L a aL^ — LpaL a p = L 2 e a pe & p can be rewritten as 


1 

4 


-L 2 e a/ 3 e A p + - ( L adl Lpp + Lp & L a p) 


(3.28) 


This is 



the “identity” or 



Thus, we can write: 


The permutation is obtained by swapping a and or 



-Zr e a/3 e^ + - (L^L^ + L^L^) , 
tapCa/i + 3 {LaaLpp + Lp&L^) . 


(3.29) 

(3.30) 


In this language, the tree-level contraction is represented as 



(3.31) 


Hence, also reinstating powers of the ’t Hooft coupling, we obtain that the term in the 
spin-chain Hamiltonian corresponding to the two-fermion —> two-fermion process is 


A 

87T 2 




(3.32) 


in agreement with the corresponding terms in (3.1). In conclusion, putting together the 
purely scalar result of Section 2, (2.20), as well as the results (3.14) and (3.32) for the 
two-fermion two-scalar and four-fermion cases, we have confirmed the complete expression 
(3.1) for the spin-chain Hamiltonian in the 577(213) sector. 


4 Conclusions 

We would like to summarise some of the key points of our paper, compare with other recent 
on-shell approaches and finally make a few suggestions for future research. 

The calculation of the dilatation operator in J\f — 4 SYM has been revisited in recent 
months using on-shell and twistor (string) inspired approaches. In [15] the complete one- 
loop dilatation operator was obtained by calculating form factors for generic single-trace 
operators using generalised unitarity, making interesting contact with earlier work of [18]. 
In particular, the integral form for the dilatation operator in [18] is mapped to a phase- 
space integral, which appears naturally in a unitarity-based approach. The calculation of 
two-loop form factors using unitarity was also employed to obtain the two-loop anomalous 
dimension of the Konishi operator in [16]. 

On the other hand, in [17], twistor-space MHV diagrams were used to find the dilatation 
operator in the 50(6) sector at one loop directly from two-point correlators, leading to the 
position-space form of the correlator as found by [9]. In a closely related approach, in [1] 
MHV diagrams in momentum space were shown to reproduce the 50(6) one-loop dilata¬ 
tion operator. In momentum space the calculation gives a single-scale two-loop integral, 
which after Fourier transform gives the expected result. Perhaps one interesting difference 
between these two MHV-based approaches is that the twistor-space computation requires 
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an additional line-splitting regularisation of the operator. Finally, in the current paper we 
have simplified the calculation of [1] considerably by directly applying generalised unitarity 
to the calculation of the two-point functions, which should have obvious generalisations to 
higher loops. 

In comparing the two main lines of approach, using form factors or the two-point cor¬ 
relators, one notices the following main points. In order to extract L-loop anomalous 
dimensions from form factors, an L-loop calculation is required, while for the two-point 
correlators in momentum space in principle 2L-loop integrals can appear. However, form 
factors also have (universal) infrared divergences which need to be disentangled from the 
UV divergences, and with increasing loop order one obtains integrals with an increasing 
number of scales. In the case of two-point correlators, one has the advantage of only 
having to consider single-scale integrals, albeit at higher-loop order in momentum space, 
and one never encounters infrared divergences. More work is clearly needed to determine 
which method is more efficient, but we think that all approaches have their own merits 
and will shed interesting new light on the problem of calculating the dilatation operator, 
and hopefully lead to a proof of the integrability conjecture. 

Let us now make some concluding comments on our hirelings and point out future 
directions. 

1. We stress that in our method no integrals are computed at any stage. We only 
identify coefficients of a single quadruple-cut integral, and from its UV divergence we read 
off the dilatation operator. Some PV reductions are performed on shell, which are also of 
algebraic nature. 

2. At one loop, quadruple cuts are sufficient to determine the dilatation operator. This 
is related to the fact that at this loop order only four fields are connected (and thus the four 
propagators we cut are always present). This circumstance is not related to the presence of 
maximal supersymmetry, and hence we can envisage obvious applications to theories with 
AT < 4 or even no supersymmetry. 

3. In the calculations presented here (as well as in [1]) we have made use of amplitudes 
with scalars and fermions. The use of gluon amplitudes remains as a future direction of 
research, and we expect these to be relevant for the study of the SL(2) sector as well as 
for single-trace operators made of held strengths in QCD [19]. 

4. Clearly the application of our method to higher loops is a crucial testing ground - 
the ultimate goal being proving integrability at higher loops (rather than assuming it). 

5. It would also be interesting to perform the SU (213) calculation with MHV diagrams, 
thus extending the approach of [1] to fermions. 

We will come back to these issues in future work. 
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